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Our study of the evolution of transmission eigenvalues, due to changes in various physical pa- 
rameters in a disordered region of arbitrary dimensions, results in a generalization of the celebrated 



DMPK equation. The evolution is shown to be governed by a single complexity parameter which 
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implies a deep level of universality of transport phenomena through a wide range of disordered re- 
gions. We also find that the interaction among eigenvalues is of many body type that has important 
consequences for the statistical behavior of transport properties. 

PACS numbers: 68.65.-k, 05.45.-a, 05.30.-d 



I. INTRODUCTION 



Recent advances in nanotechnology and quantum information theory have motivated an extensive research on the 
topic of electronic transport through disordered regions [1-3]. The random distribution of impurities in such systems 
give rise to fluctuations in transport properties from sample to sample. The fluctuations can also be observed in a 
single disordered sample under an external perturbation or a slight variation of a system parameter [4,5]. As a result, 
1 . ■ the information about the statistical behavior of transport properties is of great importance. 

A variety of transport properties can be formulated in terms of the eigenvalues of transmission matrix of the region. 
The knowledge of the statistical behavior of transmission eigenvalues is therefore very useful in the statistical analysis 
of transport properties. This motivates us to study the joint probability distribution of transmission eigenvalues. 
Previous attempts in this direction have resulted in the well-known DMPK equation which describes the statistical 
y—i | evolution of transmission eigenvalues with respect to changing length of the medium [2,6]. Various assumptions made 
■ in its derivation, however, restrict its applicability to quasi one dimensional systems or under specific scattering 
conditions [1,2,7]. This being the only analytical tool available so far to describe the full distribution of transport 
properties, a generalization of DMPK equation for higher dimensions and under generic scattering conditions is 
£^ \ required. (The other analytical method based on nonlinear sigma model provides information only about the moments 
of the transport properties) . Further the transport properties are also sensitive to changes in other system parameters 
I/"") ' besides length e.g., boundary conditions, disorder strength and dimensionality. It would also be desirable if the 
generalized equation contains explicit information about the effect of various system parameters on the distribution. 
Our study in this paper is aimed at such a generalization. 

The derivation of DMPK equation is based on the transfer matrix approach, applied to a conductor placed between 
two perfect leads of finite width. A transfer matrix relates the wave amplitudes on the right of the region to the 
left. The scattering of waves by randomly distributed impurities leads to a randomization of transferred amplitudes. 
The transfer matrix can therefore be modeled by a random matrix, that is, a matrix with all or some of its elements 
O ■ as randomly distributed [1,2,8,9]. The idea to use random matrix approach originates from the universality of the 
conductance fluctuations, observed in the metallic regime of different disordered systems. The universality suggested 
the possibility of formulating a theory which is system- independent; the suggestion motivated the use of standard 
random matrix ensembles [1,2] as models for the ensembles of transfer matrices in the metallic regime. These ensembles 
are obtained by using maximum entropy hypothesis under a single constraint of fixed eigenvalue density; their matrix 
elements are of almost same strength and the statistical behavior, being governed only by underlying symmetry, is 
universal in nature [1,2]. These ensemble can therefore serve as good models for the transfer matrix in a metallic regime 
or for quasi one dimensional conductors where the flux incident on one channel can be assumed to be transmitted with 
the same probability into all outgoing channels (known as the isotropic assumption). The DMPK equation was also 
derived by using the standard random matrix model for the transfer matrix of a small length of the disordered region. 
As expected, the equation has been very successful in predicting the behavior of transport properties in metallic 
regime or for quasi one dimensional conductors. 

The standard random matrix ensembles are not appropriate tools to model the transfer matrix of a disordered region 
under generic scattering conditions. This is confirmed by the observed failure of DMPK equation beyond metallic 
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regime, in conductors of dimensions higher than one, or under inelastic scattering conditions. The generic scattering 
conditions may cause a preferential as well as multiple interactions among channels which would lead to varying degree 
of correlations among transfer matrix elements. The equation governing the evolution of the transmission eigenvalues 
therefore should be based on a model of transfer matrices free of isotropy constraints. In this paper, we consider 
such a model and study the evolution of the transmission eigenvalues due to change of various system parameters. As 
shown later, the resulting evolution equation indicates, beyond metallic regime, the presence of eigenvalue correlations 
stronger than those suggested by the DMPK equation. This would lead to new theoretical predictions of the statistical 
behavior of transport properties. However, in the metallic regime, the bulk eigenvalue correlations suggested by our 
evolution equation are essentially the same as those given by DMPK equation. 

The paper is organized as follows. In section II, we derive the diffusion equation for the elements of the transfer 
matrix of a disordered region under generic scattering conditions. The equation is then used to obtain the statistics 
of transmission eigenvalues in section III. To maintain the flow of the discussion, only relevant steps are given in 
sections II, III; the details of the steps can be found in the appendices. The section IV deals with a derivation 
of the " Hamiltonian" formulation of the dynamics of transmission eigenvalues. The usefulness of the Hamiltonian 
representation is that it clearly reveals the hidden three body eigenvalue-interactions which are found to be absent 
in the case of DMPK equation. This is followed, in section V, by a brief discussion of the averages of the moments 
of transmission eigenvalues; the information is relevant for predictions of the distribution of the transport properties. 
The section VI contains a discussion about the scaling behavior of the transport properties. We conclude in section 
VII by summarizing our main results. 



II. STATISTICS OF TRANSFER MATRIX ELEMENTS 



We consider a disordered region of length L and width W connected to two electron reservoirs by ideal (perfectly 
conducting) leads. The scattering in the region as well as reservoirs is assumed to be elastic. The finite width in 
the transverse direction leads to the quantization of energy of the transverse part of the wavefunction. As a result, 
the scattering states at the Fermi energy satisfy the relation kp = k^ + e„ with fcf as the Fermi momentum, k n the 
longitudinal momentum {k n > 0) and e„ as the transverse quantized eigenvalue. The various k n for n = 1, 2, .., N define 
the N propagating channels. As each channel can carry two waves traveling in opposite directions, the wavefunction 
on either side of the disordered region is specified by a 27V components (corresponding to the amplitudes of N waves 
propagating to the right and N waves to the left). The normalization of the wavefunction is chosen such that it carries 
unit current. Various length-scales associated with the wavefunction divide the transport properties into three main 
regimes: (i) ballistic limit described by I > L with I as the mean free path, (ii) diffusive limit given by I < L < ( with 
£ = Nl as the localization length and (iii) insulator limit with L > £. 

The scattering properties of the disordered region are completely characterized by a 2N x 2N transfer matrix M 
with Mki = X)s=i(*) s_1 ^fe':s' the subscript s refers to the real (s — 1) and imaginary component (s = 2) of the 
element. The transfer matrix has a multiplicative property: if the disordered system is described as a sequence 
of n "scattering units" (thin slices), with transfer matrices Mi, M n , respectively, the transfer matrix of the 

disordered system is 

M (n) = M n ...M 2 M 1 (1) 
The current conservation imposes a "pseudo-unitarity" constraint on M: 

S Z M" 1 S 2 = M f (2) 

where S z is a 2N x 2iV matrix: S z = ^ ^ ^ ^ with 1 and as N x N unit matrix and null matrix, respectively. 
The presence of time-reversal symmetry in the disordered region subjects M to an additional requirement: 

SxME, = M* (3) 

with E x as a 2N x 2N matrix: T, x — 

The current conservation condition on M leads to correlations between its matrix elements. As a result, the number 
of independent elements of M is N 2 — (2N) 2 . The time-reversal symmetry along with current conservation further 
reduces the number of independent elements to iVi — N(2N + 1) [6,2]. In this section, we denote an independent 
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matrix element of M by with p as a single index, running from 1 — ► Np, and replacing three indices (k, I, s), that 
is, p = {kl; s}; the subscript (kl; s) will be used only if required for clarification. Here Np — N(2[3N + 2-/3) with 
(3 = 1 in presence of time-reversal symmetry and (3 — 2 in its absence. 

Our objective in this paper is to study the statistics of transmission eigenvalues. In principle, the eigenvalues T n 
of the transmission matrix T of a region are known as transmission eigenvalues. However, due to a simple one to one 
relationship with T n , the doubly degenerate eigenvalues A„ of the matrix B = [A + A^ 1 — 2]/4, with A = M.M\ 
are also referred as transmission eigenvalues: T n = ^ x , n = 1 — ► N [1]. The A„ can further be expressed in 
terms of 2N eigenvalues of the matrix A which exist in inverse pairs and can be denoted by x n (n = 1, 2, ...,2N) [1]: 
A„ = (x„ + x^ 1 — 2)/4. The statistics of transmission eigenvalues can then be determined by a knowledge of the 
statistics of eigenvalues x n which in turn is related to the statistics of transfer matrices. 



A. Probability density of transfer matrices and its dependence on various system parameters 

For statistical analysis of transfer matrices, we consider an ensemble of matrices M, representing a collection of 
disordered conductors of length L and defined by a differential probability dp(M) — p(M)dp(M) [2,6,10]. Here p(M) 
and dp(M) are the probability density and the invariant measure associated with the transfer matrix space (see [10] 
for the details about dp{M)). The behavior of p(M) depends on the physical properties of the disordered region. This 
can be explained as follows: The elements of M describe the overlap between various channel states on two sides of 
the region. The presence of disorder and deterministic uncertainty due to complexity of region leads to randomization 
of the elements M^. Based on complex nature of the disordered region, the randomness associated with M M can be of 
various types. Further the genetic scattering conditions can cause multiple channel interactions, resulting in varying 
degree of correlations between matrix elements. The multiplicative property (eq.(l)) of the transfer matrices also 
imposes constraints on the behavior of p(M): if pl(M) and pl {Mq) be the densities of transfer matrices M, Mq of 
a region with lengths L,Lq, respectively, the density pL'(M') of the matrix M' — MMq corresponding to the length 
L' = L + Lo should be reproducible under convolutions of pl(M) and pl {Mq), 

Pl>(M') = J p L {M'M^ l )p Lo {M )dp{M ) (4) 

The above formulation was used as a basis to derive the DMPK equation, by applying it to obtain the probability 
density for a length L' = L + SL with Lq = 5L as a small length increment. The psl(Mo) in this case was assumed to 
be a maximum entropy distribution obtained under isotropy constraints: P5l(M ) = exp(/i — z/TrM^M) with p and 
v as Lagrange Multipliers [6,10]. The assumption implies that, in DMPK case, the information about matrix-element 
correlations is contained only in the invariant measure dp(M); the latter, however, contains only those correlations 
which are required to preserve the current conservation [10]. 

The generic scattering conditions in a disordered region result in the matrix element correlations beyond those 
due to current conservation condition. We therefore need to consider a probability density p(M), free of isotropy 
constraints, with multi-channel correlations and dependent on various other system parameters besides length; (now 
onwards, we suppress the subscript length of p). In absence of any further information about the disordered region, 
the simplest and least biased hypotheses is that the system is described by the distribution p{M) that maximizes 
Shannon's information entropy 

S[p(M)} = - J p(M) lnp(M) dp(M) (5) 

under the constraints (i) p(M) is normalized, (ii) the mean < M M > and correlations < M^M^i > are fixed: 

< Mfj_ >— f dl Qa C ' < M^M^i >— [3 d g^ gC r , with C as the normalization constant and the parameters a M , b^> given 
by the system conditions. The entropy S subject to above constraints can be maximized by considering the functional 
S[p] = S[p] - {1/13) Y^n,n' < M v M ^' > -( 2 //3) J2ft a M < M v > ~ C and putting its functional derivative to zero : 



SS = - J dp(M) Sp(M) 1 + lnp(M) - I £ b^M^M^ - | ^ 
The above implies 



0. (6) 
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p(M, a, b) = Cexp 



-(1/0) b^'MnMp. - (2 /(i) ]T ciyMy 



(7) 



with C as a normalization constant and a, b as the matrices of distribution parameters a M and byy>, respectively. Note 
the symbol £^ implies a summation over distinct matrix elements only. 

The Gaussian form of p in eq.(7) results due to consideration of constraints only up to second order moments 
of the matrix elements. The availability of information about higher order moments may lead to a non-Gaussian 
behavior of p. The Gaussian assumption, however, can also be justified by the same intuitive reasoning which led to 
the Gaussian distribution of p(M) for a small length increment SL in the DMPK case [10]. The reasoning is based 
on the multiplicative property of M: Dividing the length into small segments, the transfer matrix of the region can 
be written as a sum over matrices corresponding to the length segments [11]. Assuming the matrices in the sum to 
be independent of each other, a central limit theorem for the statistics of transfer matrices was proved in [11] in the 
weak-scattering limit. For a length macroscopically small but containing many scattering units, as in the DMPK 
case, the central limit theorem justifies the assumption of independent Gaussian distributions for various My. For a 
macroscopically large length, however, the correlations among M M 's should be taken into account which intuitively 
leads us to eq.(7). 

The distribution parameters a M and byy>, being measures of the averages of the matrix elements and their cor- 
relations, are influenced by various system conditions. For example, the increase of disorder in the region leads to 
localization of waves and reduced channel-channel interaction which in turn affects the distribution of each My. Simi- 
larly the increasing dimensionality d of the region increases the total number N of existing channels, N m {kpW) d ~ l , 
as well as the probability of interacting channels. The latter increases the number of finite b parameters. The boundary 
conditions or topology of the region also affect the distribution parameters, due to their influence on the interactions 
among channels close to boundaries. As an example, consider the case of interaction between nearest-neighbors chan- 
nels only. This gives only 2N(t + 1) finite &'s, rest of them being infinite; here t, the number of nearest neighbors 
depends on the dimensionality as well as boundary conditions of the region. Further the strength of finite b parameters 
depends on the localization length Q of the region which, in turn, is quite sensitive to the dimensionality. 

Another important system parameter affecting the sets a, b, is the system length L. The dependence of a, b on L 
can be derived by using the convolution property, given by eq.(4), as a condition. For example, the L-dependence of 
6's can be seen by considering a simple case with a M length-independent. The eq.(4) gives, by using eq.(7) for p(M') 

1< 



as well as p(M) and writing b(L + SL) « b(L) + §j-SL 



„ -(1//3)(5LV -JUii M'M' ; „ f -(1//3)V ,b, llL ,[M u M u ,-M'M',] ,,,, , ,,,, , a . 

C L+SL e ^^'^^ „ »' =c L c y " J >^^y » »"p SL (M Q ) dp(M ) (8) 

where M = M'Mq 1 = M'T, Z M^T, Z . The ^f - can now be obtained by multiplying eq.(8) by H^dM'y, followed by 
an integration: 



Dct 



db 
dL 



r 2 



where / as a function of 6, f(b) = J fo(Mo)psL{M )dp(M ), and / as a function of M and b: fo(M ,b) = 
J e -(V/?) V^m^-m^m;,] Yl^dMjj. The eqs.(8,9) indicate that the variation of b with respect to L depends 
on the scattering properties of the small length increment. 

The physical properties such as current conservation and presence of time-reversal symmetry in the region also 
put restrictions on the strengths of the parameters a, b. The current conservation condition (eq.(2)) implies £J< 
MkiM* t > — < Mk,i+NM* l+N >] = 5k s Ck (with Cfc = 1 for all k < N and Ck = — 1 for k > N) which in turn connects 
various 6's (as < Mij- s 'Mki- s >= P au l ° eC )■ Similarly the time-reversal symmetry along with current conservation 
results in equality M ki = M£ +N<l+N [6] which implies b k iij; SS ' = (-l) s ~ 1 bk+Nj+N : ij;ss> = (-l) s ^bki^+N.j+N.ss' = 

{ — 1) S+S ~ 2 bk+N ,l+N ,i+N ,j+N:ss' and akl-s = ( — l) s ~ 1 ak+N,i+N:s- 

In general, the different scattering conditions can give rise to different sets of distribution parameters a, b. For 
example, let us consider the situation in three main regimes: 

(i) Ballistic Regime: The transfer matrix in this regime is almost an unity matrix. The regime can therefore be 
modeled by an ensemble of matrices M with a probability density given by eq.(7) where 

ay -> (l/2L)5yy d , , b yy, ^ {I / L)S yy> (10) 
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for all p, p! , fid- Here pd = (kk;s) with M Md = Mkk;s as a diagonal element. Note, for the case L << I, the above 
parametric-strengths correspond to an ensemble of matrices with almost non-random elements; it approaches the 
ensemble of diagonal matrices in the limit L/l — ► 0. 

(ii) Metallic Regime: For this case, (a) the overlapping between various channels is almost of the same strength, 
(b) the flow between two channels is not affected by the presence of other channels. The ensemble of matrices M can 
then be modeled by eq.(7) with distribution parameters 

a^O, , ^'-(C/iJV (11) 

for all p, p! , Hd and with I < L < £. The above set of parameters result in an isotropic density p(M) oc 
exp[— (C/2L)Tr(M.M 1 ')] which corresponds to a statistical behavior independent of system-details. The latter is 
in agreement with the observed behavior in the metallic regime. 

(in) Insulator Regime: This regime corresponds to a zero net current flow across the disordered region which implies 
an almost zero overlap between different channels states. The insulator state can therefore be modeled by the limit 

a^O, , ^-(L/C) (12) 

for almost all p,p', and, with L >> £. In large L limit, above parameters result in an ensemble of transfer matrices 
with all matrix elements going to zero. 



B. Single Parametric Evolution of p(M) 

The distribution of the transmission eigenvalues A„ can be obtained, in principle, from eq.(7) by using eq.(2.11) 
of [10] which relates M with a matrix r(A), with functions of A as its elements, and the unitary matrices U, V: 
M = Ut(X)V. An integration over matrices U, V, if possible, can then give the probability density for A. The non- 
isotropic form of the distribution p(M), however, makes the integration route very complicated. This motivates us to 
seek an alternative route. As discussed below, we reduce the technical complications by formulating the non-isotropic 
problem in the same form as that of the isotropic problem. This also helps in identifying a single parameter which 
governs the evolution of multi-parametric p(M). 

We proceed as follows. A perturbation of the disordered region due to a change in impurity structure or other 
system parameters perturbs the matrix elements M M and, consequently, the probability density p(M,a,b). Due to 
its Gaussian form, a change in p due to variation of M M can be well-mimicked by the change due to the distribution 
parameters a, b (see appendix A for the derivation): 



Lp = Tp 



(13) 



where L and T are the operators in M-space and parametric space, respectively, 



L = E 



d 



d 



2 8M U 



da,, 



(14) 
(15) 



with 



(16) 



Tfifi' 



G 



fifj,' 



7*w - y~] W" 



(17) 



Here = 1 + S ik Sj l if p, = (kl; s), p! = (ij; s'), and, S xy = 1 if x = y and S xy — for x ^ y. The normalization 



constant C is chosen such that TC = X^u(7 + (^P 



bfj,fj,))C (see appendix A for details). The parameter 7 



is arbitrary and marks the end of the transition. Note, in eq.(13), the derivatives of p with respect to different 
matrix elements are independent of each other but the parametric derivatives are correlated. By considering the 
above combination of derivatives, therefore, we transfer the information regarding the correlations between elements 
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of matrix M to the elements of matrix b. This helps in reducing the M-space operator L (eq.(14)), in the same 
form as the one in the isotropic case; as mentioned above, this reduction is useful for technical reasons as well as for 
comparisons with earlier studies. 

Contrary to the single parametric evolution of p in the isotropic case, the evolution of p in eq.(13) is governed by 
a large number of parameters. However, as discussed in appendix B, seemingly multi-parametric diffusion can be 
reduced to a single parametric evolution by considering a transformation (a, b) — ► y, with y as a set of parameters 
yj, j = 1, 2, .., N, which maps the parametric space operator T (cq.(15)) as T[a, b] — ► T[y(a, b)] = Here TV is the 
number of (a, b) parameters with finite strengths. As a result, eq.(13) can be written as 



E 



d 



Q p 

1 7^ p 



2 8M fl 



dp_ 

dY 



(18) 



where Y = y\. The eq.(18) is in a form of a standard Fokker-Planck equation with "particles" undergoing a 
Brownian motion in "time" Y. The evolution approaches a steady state in limit — > or Y — > oo which occurs 
when ffnfufj,' — * or, equivalently, — > 0, & w < — > 7^^' and, therefore, p — > e -(7/ 2 )' rrAfM+ (from eq.(7)). The latter 
is also the solution of cq.(18) in limit — > oo; the steady state limits of eq.(7) and (18) are, therefore, consistent 
with each other. 

The parameter Y, appearing in eq.(18), is a function of various parameters and b^ which are governed by 
underlying complexity of the system; Y can therefore be termed as the complexity parameter, 



Y = ^2 / da M z n x + E / •W X + constant 



(19) 



where summation is implied over the distribution parameters with finite values only, and, X — z^fn + 
J2fi n' z nn' fvit'] -1 ■ Here the functions z^, are arbitrarily chosen such that the ratio 



H n,n> 



X (20) 



is a complete differential; the details of the eqs. (19,20) are given in appendix B. 

It is worth emphasizing here that the system information in eq.(18) enters only through the parameter Y. The 
distribution parameters a, b being system-dependent, Y in turn is a function of the system parameters e.g. length, 
disorder, dimensionality and boundary conditions etc. A variation of any one of the system parameters e.g. length L 
can therefore change Y but, note, in general Y ^ L. (In other words, even if L is the only system parameter subjected 
to change, the parameter governing the evolution is Y, a function of other system parameters besides L). However 
the case Y cx L can occur if the distribution parameters are assumed to depend only on the length of the system 
(ignoring the dependence on other system parameters). For example, consider the case with 



for all /z, p! with q as a constant. Using the values in eqs. (16, 17) gives = 0, / w = 2(qL)~ 2 (-~fqL — 1) and therefore 
X^ 1 = 2N p(qL)~ 2 [^qL — 1]; here we have chosen = 1 as that makes the ratio (20) a complete differential. Note 
the variation of length L is this case changes only Np parameters b^* , leaving all other parameters unchanged. This, 
alongwith substitution of X, in eq.(19) gives 



Y 



' J 2(7^-1) = ~^ l0g|1 ~ 7qLl + COnStant ' (21) 



As obvious, the system length L ~ (79) 1 gives Y — > 00 and therefore steady state of the evolution: p — > 

e (-7/2)TrM.M t 

The dimensional sensitivity of Y can also be explained by a simple example. Let us consider a disordered region 
with interactions between nearest-neighbor channels only, that is, by choosing 

a (U = ^ w , -^rc^S^, (22) 

where = 1 if M M corresponds to an interaction between same channels (i.e. for p = pd = (kk;s)) or nearest 
neighbors k,l (p = (kl; s)), c M = for all other p values, (thus ^ M c M = 2N(t + 1) with t as the number of nearest 
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neighbors). As channel-channel interaction depends on the localization length £, it is reasonable to consider r = Q/L. 
For the above a, 6-values, eqs.(16,17) give = (7 - 2r)<5 wid , / w = 2rc M (7 - rc M ). Choosing again = z^ = 1, we 
get X^ 1 = 2N[j + 2-yr - 2(t + l)r 2 . Substituting the above in eq.(19) gives 

Y = 2(iN j drX 



I3N , 2{t+l)C,+q+L 

" ln ^7 7\7 — ^ + constant (23) 



27(t + l) 2(t + l)C-g_£ 



with <?± = y 7 7 2 t 2 + (t + 1)7 ± -ft. The parameter being a function of dimensionality as well as boundary conditions 
of the region, its presence in eq.(23) implies the dependence of Y on the conditions. However the most significant 
dimension-dependence of Y comes through the localization length £ which varies with dimension of the region (see 
[12]). Further, as steady state of the evolution (Y — > 00) occurs at length scales Q ~ 7L, the approach to equilibrium 
is sensitive to system-dimension. 

The parametric space transformation (a, b) — > y maps the probability density p(M, a, b) to p(M, y(a, b)). As a result, 
p depends on various parameters yj , j = 1 — > N with JV as the number of non-zero elements in set (a,b) . However the 
diffusion, generated by the operator L in the matrix-space M, is governed by Y = y\ only; the rest of them, namely, 
yj, j > 1 remain constant during the evolution. Note it is always possible to define a transformation from the set 
(a, b) — > y with yj, j > 1 as the constants of dynamics generated by L. This can be explained as follows. A matrix 
element, say M M , describes how a basis state ipi interacts with state ipj through M. This results in dependence of the 
matrix element correlations and, thereby, of the set (a, b), on the basis parameters e.g. basis indices. As the basis is 
kept fixed during the evolution, the suitable functions of basis parameters can be chosen to play the role of yj, j > 1. 
(Note a similar transformation has been used to obtain a single parametric evolution of multi-parametric Gaussian 
ensembles of Hcrmitian matrices; see [13] for details). 

The eq.(18) forms the basis of our results obtained in next few sections. As clear, it is based on the Gaussian 
assumption for the probability density of the transfer matrices. The complexity of the region may however give rise 
to non- negligible higher order matrix element correlations which leads to non-Gaussian behavior of the p(M). In this 
context, it is useful to note that the evolution of the density p(H) of the Hermitian matrices H is also known to 
be governed by an equation similar to eq.(18) (with M replaced by H and a different form of Y); the p(H) can be 
Gaussian or non-Gaussian (see [18] for the proof). Following similar steps as used for the non-Gaussian p(H), the 
eq.(18) can also be proved for a non-Gaussian density of transfer matrices e.g. for p(M) oc e^ M \ with f{M) now 
containing products of r matrix elements, r = 1 — > n. 



III. EVOLUTION OF TRANSMISSION EIGENVALUES DUE TO CHANGING COMPLEXITY 

The single parametric distribution P(X, Y) of the transmission eigenvalues can, in principle, be obtained by sub- 
stituting M — Ut(X)V in the solution p(M,Y) of eq.(18), followed by an integration over matrices U, V. However, 
again, as in the case of eq.(7), the integration route offers many technical difficulties. This again motivates us to use 
the differentiation option which has, as discussed later, other advantages too. 

In this paper, we use the evolution of M, given by eq.(18), to show a single parametric diffusion of the transmission 
eigenvalues. The steps are given as follows. The eq.(18) gives the first and second moments of the matrix element 
components Mki :s = as (appendix C) 

SM klls = -iM khs SY, 5M khs 5M mn;s , = [36 km 6i n S ss ,SY (24) 

with / implying an average of / over random noise. (For clarification, we again refer a matrix element component as 
Mfcj ;s ). The above moments can further be used to determine the moments of matrix elements A mn — J2k=i M m kM* k 
(see appendix D for details), 

5A mn = —2-fA mn SY 

5A mn 5A* mn = 2 [/? + (2 - f3)S mn + (4 - f3)S mn+N )] (A mm + A nn )8Y (25) 

with (3 — 1 in presence of the time-reversal symmetry and (3 = 2 in its absence. 

As mentioned before, the matrix A has 2N eigenvalues x n , 1 < n < 2N which form inverse pairs. For clarification, 
let us label the eigenvalues such that x n+ N = x^ 1 . Using 2 nd order perturbation theory for Hermitian matrices, the 
change Sx n of the eigenvalues x n due to a small perturbation 5 A of the matrix A can be given as 



7 



Sx n = 6A nn + J2 + 0((5A mn f) 



The above equation can then be used to obtain the first and second moments of the eigenvalues x n , n = 1 — > N 



5x n = F n (x)5Y, Sx n Sx m = 8x n S nm 6Y 



(26) 



with F n = -2-fx n + 2J2 m ^ n [P + (4 - (3)S m ^ n+N ](x n + x m )(x n - x m ) 1 (appendix E). The cq.(26) along with the 

relations A„ = , 5X n = ^f^-5x n gives the first and second moments of the transmission eigenvalues A„ 

(appendix E), 



5\ n 5\ m = 8 A " (A " + l) S mn SY, 5X n = —AE n (X n )SY 



(27) 



where 



E n (X n ) = —[^X n y/X n (l + Xn)- (2A„ + 1) - /3 V 



A n (l + A ra ) 

\ — \ * 



(28) 



with x n = 1 + 2A„ + 2yJ\ n (\ n + 1). The higher order moments of 5X n vanish to first order in 5Y. 

The information about moments of the eigenvalues A n can now be used to obtain their evolution equation. The 
theory of Brownian motion [14] informs us that the joint probability distribution P({A}) for the eigenvalues A„ evolves 
with increasing Y according to Fokkcr-Planck equation, 



op i A d 

n=l 



l d(SX n 5Xm)P , , jr-. p 

2 dX n + (- 5A ») P 



N a 
^dX n 



n=l 



a(x- 1 A„(l + A»)P) 
9A„ 



(29) 



Equation (29) describes the evolution of the transmission eigenvalues for a disordered region with respect to its 
changing complexity [15]. The reasons for a change in complexity can be manifold, for example, due to change in 
length or disorder of the region, the scattering and boundary conditions. As a result, eq.(29) is not only valid for 
higher dimensions and beyond weak scattering limits, it can also describe the influence of various system parameters 
on transport; the information about all of them enters in the probability distribution through a single parameter Y. 
Note Y can also be expressed as a function of the localization length ( measured in units of system length L. This 
can be explained as follows. A single parametric formulation of the diffusion in the matrix space M implies a same 
parametric dependence of the evolutions in the eigenvalue as well as eigenfunction space. This allows a complexity 
parameter formulation of the eigenfunction correlations e.g < V'nMV'n^') >= f{r, r'; Y) (with ip n (r) as n th eigenstate 
of the transfer matrix at space point r of the disordered region). As the localization length ( can be expressed in 
terms of eigenfunction correlation f(r, r'Y), this results in Y dependence of (. 

It is worth comparing eq.(29) with DMPK equation. The latter describes the evolution of transmission eigenvalues 
due to changing length of the region in the diffusive regime and for quasi one-dimensional systems. The DMPK 
equation has been obtained by assuming the elements of the transfer matrix of a small length L of the region as 
independently distributed: pdmpk(M) oc c -("'/2i)TrMM t with a = /J_/v + 2 - /3. Note that pdmpk is a special limit of 
eq.(7), equivalent to a choice of all a M — ► 0, b^^< — > (al/L)S^^. Consequently the complexity parameter Y for pdmpk 
turns out to be same as eq.(21) with q — (al)" 1 . In the limit L < Nl, Y reduces in a form similar to the evolution 
parameter of DMPK equation: Y — ► L/2al. However the right side of eq.(29) is still slightly different from that of 
the DMPK equation. By further considering the strongly transmitting channels limit, namely, A„ << 1 or x n w 1, 
eq.(29) reduces to DMPK equation 



N 



dP d 

= 4 E 



dY 



8X n 



^(l + X n )P]__ {2Xn + l)p _^X n (l + X n )P 



8X n 



X n — X r: 



(30) 



Thus, for strongly transmitting channels and with length as the only changing parameter, eq.(29) describes the same 
statistical properties of the transmission eigenvalues as the DMPK equation. However, as eq.(30) indicates, if the 
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system parameters other than length are subjected to change, eq.(29) has a different evolution parameter even in 
strongly transmitting channels limit although its A-dependcnt terms remain same as in the DMPK equation. The 
latter similarity is useful for the following reason. As many results for the length-dependence of transport properties 
e.g. conductance using DMPK equation have already been derived, a replacement of length by Y in the known results 
can then give us the dependence of properties on other system parameters too (in limit A„ << 1 only). 



IV. HAMILTONIAN FORMULATION OF THE EIGENVALUE-DYNAMICS 



The distribution of transmission eigenvalues and their correlations can be used to determine the statistical behavior 
of various transport properties. For example, the conductance G is given by G — 5Z n (l + A„) _1 in the Landauer 
formulation (in units of dimensionless conductance Go = 2e 2 /h) [1]; consequently, its distribution is related to that 
of the transmission eigenvalues: 



AT N 

P G (G) = S(G-J2 TTT-) P ({ A }) II dA 

J n=l + A ™ j=l 



(31) 



A knowledge of the solution of eq.(29) is therefore very desirable. Although the solution is not known so far, still 
we can gain a deeper insight about eigenvalue correlations by writing eq.(29) in a Schrodinger equation form. The 
required steps are as follows. A transformation of the eigenvalues A„ = sinh 2 (log/i„) reduces eq.(29) to a Fokker-Planck 
equation with constant diffusion coefficient, 



dP 



N 



d 



i)Y S a 

n=l 



dp | 0nfc n ) 

d/J, n dfi n 



(32) 



where Q(/i„) = - J2 n >m lo §l A ( r ™> r ™) - Z 3 " 1 En[ lo g( sinn2r ™) _ 7M™sinhr„ +r n ] + (N(N - l)ln2)/2 with A(r„, r m ) = 
(cosh2r„ — cosh2r m ) and r n = log/x„. By using P({[i}, Y) = e~ /3n / 2 ip({ii}, Y), eq.(32) can be mapped to a Schrodinger 
equation 



&Y 



= Hip 



(33) 



with H as the " Hamiltonian" governing the eigenvalue dynamics: H({/j,}) — J2 n + ^({m}) where 

2 



V(W) 



■IE 



( dQ 
2 \dfi n 



d 2 n 



(34) 



=y- 



\{^l^ 2 n ~ l 2 ^i - 3) + 7^4 coth2r„ + . \ 
4 sinh 2r„ 



/3(4 + 7 ^#^-W-2)- Sil!i " 2/ ' 



A(r„,r m ) 



A 2 (r„,r m ) 



j,jjtm,n 



sinh 2r r , 



'^n •> I'm ) ^ {j"n ; Tj ) 



(35) 



The Schrodinger equation formulation of eq.(29) reveals an important aspect of the eigenvalues, namely, the presence 
of a three body interaction terms; (note it could have been converted into a two body term in absence of the prefactor 
/i 2 ). It is well-known that the "Hamiltonian" appearing in the Schrodinger equation formulation of the DMPK 
equation contains only two body interaction terms [1,16]. Thus it appears that the presence of significant multiple 
channel interactions leads to three body eigenvalue correlations while the isotropic channel limit (used in DMPK case) 
restricts it to pairs of eigenvalues. The different "Hamiltonian" in the two cases result in a significant difference in 
their states ip and therefore P({/x}). As a result, the theoretical prediction of the physical properties under multiple 
channel interactions are expected to differ significantly from the predictions based on isotropic channel considerations; 
the suggestion is in accordance with the already known failure of DMPK equation beyond metallic regime or for two 
and higher dimensions. 
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V. EVOLUTION OF AVERAGES 



The presence of many body interactions among eigenvalues makes it difficult to calculate P({\}) under general 
scattering conditions. In the metallic regime, however, the fluctuations of transport properties are much weaker than 
the average behavior and a knowledge of their averages (instead of full distribution) is sufficient. The evolution of 
transmission eigenvalues can then be used to study the effect of changing complexity on the average behavior of their 
various functions. Consider a function i^({A}) not explicitly dependent on the parameter Y. The multiplication of F 
on both sides of eq.(29), followed by partial integration, leads to following behavior of < F >= / F({A})P({A})d{A} 
with < . > implying an ensemble average: 




d F dF 



(36) 



For example, evolution of the moment Tf — E n (l + A„) 1 ] p , due to changing complexity in the strongly trans- 
mitting channel limit, can be obtained by substituting F = Tf in eq.(36) and using the limit A„ << 1 : 



9<T *>=-2p 



13 < Tf +1 > -(/? - 2) < Tf _1 T 2 > -2(p - 1) < (T 2 - T 3 )Tf~ 2 >1 (37) 



dY 

where notation Tjf implies y th power of the moment T x = 5Z n (l + \ n )~ x ■ A substitution of p = 1, 2 in eq.(37) then 
gives the complexity dependence of the average conductance < G >~< T\ > and its variance a 2 —< T 2 > — < Ti > 2 , 
respectively, 



d < G > 
dY 

~dY 



= -2 [0 < G 2 > -(0-2) < T 2 >] (38) 
= -4 [/?(< G 3 > - < G >< G 2 >)-(/?- 2)(< M 2 >< G > + < M 2 G >) - 2 < (T 2 - T 3 ) >] (39) 



In general, the studies of fluctuations in transport properties about their average behavior require a knowledge of 
the averages of various combinations of moments e.g. < T P T% >. A substitution of F = T^T^ in eq.(36) leads to a 
hierarchy of coupled set of equations for the average behavior of the moments. The details of these calculations with 
DMPK equation as a basis are given in appendix B of [10]; the replacement of the length parameter by complexity 
parameter in eq.(B8) of [10], gives us the evolution of < T P T^ > with respect to changing complexity of a system in 
strongly transmitting channel limit. 

The coupled form of eqs. (38,39) or eq.(B8) of [10] makes it difficult to obtain the exact behavior of the moments, 
even in A„ << 1 limit . However, by assuming the number TV of the channels very large, the equations can be solved 
by using TV -1 expansion of the moments, see [10] for details. For example, the evolution of Tf reduces to, in the 
leading order of N, 

= -20p (Tf +1 ) + o(Nn (40) 

In ballistic limit, the eigenvalues A„ — > (for n = 1 — > TV) which gives < Tf >— N p . For ballistic case as the initial 
state of the disordered region, that is, lim Y — > Y < Tf >= N p , the solution of eq.(40) can be given as 

< Tf >= N p (l + A)- p + o(N p - 1 ) (41) 

with A = 2/3TV|y — Y \. Following cq.(21), A = qL/2 for almost same strength of interactions between any two 
channels and length as the only changing parameter (for jqL < 1). The behavior of < Tf > with respect to the above 
A coincides with one obtained by using DMPK equation for which q^ 1 = (j3N + (3 — 2)1 (see eq.(3.31) of [1]). The 
substitution of p = 1 in eq.(41) gives Y dependence of the average conductance < G >=< T\ >« TVA -1 + o(TV°) 
(in large TV limit). The average conductance of a wire therefore decreases linearly with increasing length which is in 
agreement with the already known results [1]. 

Equations (22,23) can be used to study the effect of changing complexity on the conductance fluctuations in the 
diffusive regime. For example, the weak localization correction SG =< G > —(1 + A) -1 to the average conductance 
as well as the variance a 2 =< G 2 > — < G > 2 can be obtained by expanding the moments to order TV -1 : 

6G = (1 - 2//3)A 3 /3(l + A) 3 + o(/V- 1 ) (42) 
Var G = (2/15/3)[l - (1 + 6A)/(1 + A) 6 ] + o^ 1 ). (43) 
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(Above results were calculated in [10] for length variation. Due to similarity of eq.(30) with DMPK equation in 
limit A„ << 1, same results can be used to describe the complexity dependence of conductance fluctuations just by 
replacing L/l by A). The diffusive limit A — > oo of the above results give SG = (1 — 2//3)/3 and Var G — (2/15/3) 
which is in agreement with expected universal behavior of the conductance fluctuations in this limit (see [1] for the 
latter). 

As mentioned above, eqs. (37-43) are valid only in the limit A„ << 1. A detailed analysis of the dependence of 
transport properties on A, beyond strongly transmitting channel limit, is still under investigation; the results will be 
given elsewhere. 



VI. SINGLE PARAMETER SCALING OF TRANSPORT PROPERTIES 

A distribution P(X) of a physical quantity X that depends on system size N and a set of system parameters {alpha,/-} 
obeys one parameter scaling if it is a function of only X and one scale dependent parameter: P(X; L; {ctfc}) = f(X; Y) 
[21]. Thus the system-information appears in the distribution only through the scaling parameter Y = Y(N, {«&})■ 

Following the discussion given in sections II, III, it is clear that the distribution of transmission eigenvalues satisfies 
the above requirement, that is, P({A„}; N; a, b) = F({A„}; Y), with complexity parameter Y = Y(N, a, b) playing the 
role of scaling parameter. As a consequence, the distribution Pf{F) of any transport property F({A}), would also 
follow a single parameter scaling: P F (F; Y) = J 5{F - F({\}))P({\}; Y) ]jf =1 d\j 

The Y governed evolution of the distribution Pf{F; Y) can be obtained with the help of eq.(29). The evolution as 
a function of Y, however, is abrupt in large AT-limit (see, for example, eqs. (40)); a smooth crossover can be seen only 
in terms of a rescaled parameter A where 

A = ^ (44) 

with D as the local mean level spacing of transmission eigenvalues A„. (Note the relation A = 2(3N\Y — Y \, given in 
section V, is valid only in the strongly transmitting channel regime). 

The parameter A governs the statistical behavior of the transport properties. As discussed in section I, the 
dimensional-sensitivity of the interaction between various channels can influence both Y and D. This in turn leads 
to dimensional-dependence of A. The statistical behavior of transport properties e.g. conductance is therefore ex- 
pected to be different for higher dimensions. The size-dependence of transport properties and influence of scattering 
conditions on them can also be explained by Y and therefore A-dependence on these system conditions. 

The one-parameter scaling behavior of the distribution P({A}) implies the existence of a universal distribution 

P*{{^\) = jy^oo A) at a critical point which is fixed by the critical value A* = A(N). Thus the 

size-dependence of A plays a crucial role in locating the critical point of statistics. It is possible that various system- 
conditions in a disordered region may result in different A^-dependence of \Y — Y"o and D, say, \Y — Yq\ oc N a and 
D oc N v which gives A oc N a ~ v . In regions connected with leads of finite width (implying finite N) , a variation 
of size N, therefore, leads to a smooth crossover of statistics between an initial state (A — > 0) and the equilibrium 
(A — > oo); the intermediate statistics belongs to an infinite family of ensembles, parameterized by A and given by the 
solution of eq.(29). However, for system- conditions leading to a = rj, the statistics becomes universal for all sizes, A 
being A-indepcndent; the corresponding system conditions can then be referred as the critical conditions. 

As clear from above, a variation of system conditions (i.e. changing a, 77) for a given finite system size, would also 
lead to a crossover of statistics. In large A^ limit, however, the variation causes a transition of statistics. Based on 
relative values of a and 77, the statistics for large system sizes can be divided into three different types 

Case a < r\: This corresponds to A — > in limit N — > 00. The system under variation of such conditions rapidly 
approaches its initial state. 

Case a > 77: This gives A — > 00 for infinite system sizes; as a result, even a small perturbation of the system will 
cause an abrupt transition from the initial state to the equilibrium. 

Case a = f3: Here A being size-independent, the case corresponds to the critical point of evolution. Further, due to 
finite and non-zero A, the eigenvalue statistics for this case is different, even in limit N — > 00, from both, the initial 
state as well as the equilibrium. Note, however, the statistics is same as the one at the critical point of finite systems. 

As clear from preceding sections (III-V), the A- formulation of transport properties (e.g. eqs. (29-43)) is based on 
eq.(18). Although this paper contains the derivation of eq.(18) only for Gaussian p(M), but, for the reasons mentioned 
in section II, we expect it to be valid for non-Gaussian cases too. This suggests the application of our results for a 
wide range of disordered regions. Also note that eq.(18) and eq.(29) are valid under inelastic scattering conditions too. 
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However the connection between the transport properties and A's needs to be revised. For example, the Landauer's 
formula G = 5Z„(1 + A„) _1 is not applicable in the inelastic regime; the results for Conductance fluctuations given 
in section V are therefore not valid for that regime. 

The validity of eqs. (18,29) for non-Gaussian cases as well as inelastic conditions suggests the existence of a A- 
formulation of transport properties for regions with generic scattering and disorder conditions. Although the functional 
form of A is expected to be different under generic conditions, the critical statistics would still occur at system 
conditions satisfying the critical criteria a = (3, thus leading to an A^-independent A. It should be stressed, however, 
that the critical criteria may not be fulfilled for all disordered systems; such systems will not show any transition in 
the statistics of transmission eigenvalues. 



In summary, we find that, under generic scattering conditions, the eigenvalues are subjected to many body interac- 
tions besides logarithmic repulsion. As a consequence, the correlations among eigenvalues and thereby the fluctuations 
of transport properties beyond metallic regime are significantly different from those within the regime. An explicit 
formulation of the fluctuations can be obtained by the knowledge of the solution of eq.(29) (or eq.(33)). Although a 
complete solution is not known so far, nonetheless eq.(29) reveals, for the first time, a very important characteristic of 
transport phenomena through disordered regions, namely, their dependence on a single complexity parameter Y . This 
implies that the disordered regions with different system parameters e.g., boundary conditions, scattering conditions 
etc. will show same fluctuations of transport properties if their parameters A correspond to a same value. This 
implies a deeper level of universality lying underneath the transport phenomena, even beyond metallic regime (where 
it was known to exist so far). A knowledge of universality is useful from a practical viewpoint, besides other reasons. 
This is because the information about the transport properties of a simple disordered region (obtained by some other 
technique) can now be used to determine those for a complicated region, just by comparing their complexity parame- 
ters. Further the dependence of distribution of the transmission eigenvalues on various system parameters through a 
single parameter Y also indicates the validity of one parameter scaling theory of localization for statistical behavior 
of the conductance (also for those transport properties which can be formulated as the functions of transmission 
eigenvalues) . 



VII. CONCLUSION 
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APPENDIX A: PROOF OF EQUATION (13) 



Using eq.(7), a derivative of po — p/C with respect to = Mki- S can be written as 




(Al) 



where b 




Further 



where = 1 + 5 i3 Ski if P = {kl; s}, p! = {ij; s'}. The eqs.(Al,A2,A3) can be combined to give 




(A2) 



(A3) 




(A4) 



Similarly 
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d 2 Po 

dM, 



M 



, y | ) {W^al - b^) P0 + ^J2 b ^' M »' ( 2fl M + E W^M" I 



/? 2 



m' 



(I) (2/r _ 6wt)po _ I ( S <w w w + 2 E a ^ 

\fi',fi" M M M' 



dpo 



Using above equalities, 



E 

M 

= E 



<9M U 



P_dpo_ 

2 9M U 



+ 7^m Po 



iG^b^-^p- - E g m'm" 5 mm' 6 



dp 



E« 



9a,, 



mm 



m' 



9a,,/ 



C/3A) 



(A5) 
(A6) 

(A7) 
(A8) 



where Cyj = [7 + ( 2 /^ — ^w*)] ■ Now interchanging the indices in the terms inside the square brackets in 
eq.(AlO) gives 



E 



dM„, 



(5 dp 



2 dM„ 



^2fw'au , +E^ 

mm 



9po 
<9a„ 



Cppo 



where / M and / w ' are given by eqs.(16,17). Multiplying eq.(A9) by C and using Cpo — p, we get 



E 



9 



9M U 



/3 9 P 



2 9M, 



5> 

ft, ^ 



db 



dp dC 
P- 



db 



'mm' 



+E/m 



dp dC 
da^ 9a M 



(A9) 



(A10) 



The eq.(13) can now be obtained by choosing the normalization condition on p such that C satisfies the relation 
/mm' db~~7 + S M /m9^ = G /3 G 



APPENDIX B: SINGLE PARAMETRIC FORM OF T: DERIVATION OF EQ.(18) 

Consider a transformation from the TV-dimensional {a, 6}-space to another parametric space, say y-space, consisting 
of variables i/j, i = 1 — ► N which reduces T(a, b) (cq.(15)) into 



T(y[a,b])p= ^\ y2 ,.., yM 



(Bl) 



. The conditions to determine desired transformation can be obtained as follows, By using partial differentiation, T(b) 
given by eq.(15) can be transformed in terms of the derivatives with respect to y, 



where 



7> = E^ 



ife = E^l!r + E/< 



dp 
dyk 



da b 



dyk 
db^> 



(B2) 



(B3) 



M ' M,M' 

The eq.(B2) can be reduced in the desired form of eq.(Bl), if the transformation b — > y satisfies following condition: 

t k =6 kl for k = 1 -> N (B4) 

According to the theory of partial differential equations (PDE) [20] , the general solution of a linear PDE 

N 



v ---> dZ 

'Y^P i {x 1 ,X 2 , ..,Xjf) — = R 



dxi 



(B5) 
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is F(m,U2, -;U n ) = where F is cin arbitrary function and %2i 

Z) = Cj (a constant), i = 1,2,.., n are 

independent solutions of the following equation 



dx\ _ dx 2 _ dx k _ dx^ _ dZ 



(B6) 



Pi P2 Pk Pn R 

Note the function F being arbitrary, it can also be chosen as 

F = Y,{u j -c j ) = Q (B7) 

3 

The equations for various yj in the set of eq.(B4) are of the same form as eq.(B5) and, therefore, can be solved 
as described above. Let us first consider the equation for y\\ its general solution can be given by a relation 
F(ui, U2, .., U/y) = where function F is arbitrary and Uj arc the functions of iV parameters of set {a, b} such 
that Uj({a,b},yi) = Cj (with Cj's as constants). The functions Uj are the independent solutions of the equation 

db uu i da u , , , 

-7^- = .... = ~r = = dyi (B8) 

where the equality between ratios is implied for all possible indices /i, u' (with the total number of ratios as N + 1). It 

is easy to see that each of the above ratios is equal to ^ — - — where z^, z^t are arbitrary functions. 

The eq.(B8) can therefore be rewritten as 

dyi = "y — z >f ,+v z f (B9) 

A solution, say ui of eq.(B9), or alternatively eq.(B8), can now be obtained by choosing the functions in the set z 
such that the right side of the above equation becomes an exact differential: 

Ml = yi - J da^z^ X - J db^ z^ X = constant (BIO ) 

A 1 M.A 1 ' 



where X = + , z^'f^'] -1 . The general solution for yi can therefore be given by a combination of all 

possible functions u obtained by using arbitrary set of z-functions. It can be shown that each such solution differs 
from the other only by a constant: Uj = ui + constant; (this is due to equality of the two ratios obtained by choosing 
two different sets z of the functions). The y\ can therefore be written as follows, 

Vi =^2 J da M z n x + X! J d6 ^' X + con stant (Bll) 

H fj,,fj,' 

which gives eq.(19) for Y = y\. 

The set of equations (B4) can similarly be solved for other yj (j > 2). For example, the solution for y k can be given 
by the function Fk(v\, ..%) = where Vj({a, 6}, y k ) = constant are the independent solutions of following equality 

db^ _ da M _ dy fc 



A solution, say vi, of eq.(B12) can now be given as 



(B12) 



da^z^ — constant (B13) 
where the set z is a set of arbitrarily chosen functions which satisfy the condition 

As obvious, one possible choice for z functions satisfying the above condition is z^i = 0, z M = for all a, which 
gives y k — constant. As each solution of eq.(B12) is different from the other only by a constant, the general solution 
for yk, k > 1, can now be given as 

Vk = ^2 dfy^'-W +^ / da M z M + constant. (B15) 
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APPENDIX C: PROOF OF EQUATION (7) 



The eq.(7) can be obtained from eq.(18), by following standard routes of calculation of the moments from a Fokker- 
Planck equation; see [17,14] for details. For completeness sake, we describe here one such route. 

Let p(M — SM; Y) be the joint probability density that the elements of transfer matrix will be at the positions 
M — SM = {M^ — SM^} at parameter value Y. At Y + 5Y, let the positions of the elements change to M = {M M }. 
The changes SM^, for all p values, are random variables and different, in general, for every member of the ensemble. 
The probability density p(M; Y + 5Y) can be expressed in terms of p(M — SM; Y): 

p(M; Y + SY) = J d(5M)tp(M - SM, Y\M, Y + 5Y)p{M - SM; Y) (CI) 

where tp is the probability that the position of the transfer matrix will change from M — SM to M in a parametric 
interval SY. Expanding p on both sides of eq.(Cl) in a power series of SM and SY around p(M, Y), one gets 



p{M,Y) + ^5Y + o{{5Yf) 



d(SM)i> 



1 3^ 

———(SM^SMij.s,) + o((5M) 



dMp 2 Jr> , dM„M ij;s , 

/J ^ khj;ss' ^ J ■ 



(C2) 



Now by using the definition (SM) n = J (SM) n ^(M^SM, Y\M, Y+5Y)A{5M) and the equality / ip(M-8M, Y\M, Y+ 
SY)d(SM) = 1, the eq.(C2) can be reduced to 



^.^[in^Jg^-SPHa^ (OS) 

By comparing eq.(C3) with eq.(18) one obtains eq.(24). The latter can also be derived by a direct integration of 
eq.(18); see [17] for details. 



APPENDIX D: PROOF OF EQUATION (25) 

The matrix A is defined by A = M.M with its elements A mn — Ylt=i M m kM* k . The average of a small change 
SA over random noise can then be expressed in terms of SM: 

2N 



5A mn = [M mk SM* k + M* nk JM^\ (Dl) 
fe=i 

27V 

= -7 E i M ^ M nk + M* nk M mk ] = -2 7 A mn 5Y. (D2) 



k=i 



Here eq.(D2) is obtained from eq.(Dl) by using the relation (24). 

The current conservation condition on the transfer matrix M, given by eq.(2), introduces various relations between 
its matrix elements: 



N 



J2 WnkM* mk - M nk+N M^ k+N ] = gS nm (D3) 
fe=i 

N 

[MnkM* n+Nk - M nk+N M* n+Nk+N ] = (D4) 



fe=i 



with g = 1 if n < m and g = — 1 for n > m. The presence of time- reversal symmetry alongwith current conservation, 
further implies that M mk — ^* m+N \\ k+N \ (M = x if x < 2N, \x\ = x — 2N if x > 2N). As a result, the second moment 
of SA mn depends on the presence or absence of time-reversal symmetry in the disordered region and is different for 
the cases m — n, m ^ n and m = n + N. 
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1. Case ra ^ n, m / n + N 

2N 

SA mn 6A* mn - Yl i M ^SM* k + Mp3 [M^SM nl + M nl SM^} (D5) 



k,i=i 

2N 



£ \M mk M* ml {5M* k 5M nl ) + M: k M nl {8M mk 5M* ml ) 



k.i=i 



+ M mk M nl {5M* k 5M* ml ) + M* nk M* ml {5M mk 5M nl )\ (D6) 

Following eq.(24), we have 

(SM jk SM rl ) = i Sl+S2 - 2 (5M jk , Sl 8M rkS2 ) = 

si,s 2 

(<5M* fc 5M r( ) = ^ (-lr^^+^I^M^M^J = 20S jr S u SY (D7) 

si,s 2 

Consequently the last two terms inside the square bracket [] in eq.(D6) do not contribute in this case and, with help 
of eq.(24), we get 

27V 

5A mn SA* nn = 2f35Y ]T [5 M M mk M* ml + 5 kl M* nk M nl ] (D8) 

k,l 

= 2l3[A mm + A nn ]5Y (D9) 

2. Case m = n + N 
Using current conservation condition given by eq.(D4), we can write 

N 

+Nk 

(D10) 

k=l 

The above equation gives 

N 

5A rwn+N 5A* n n+N = 4 [M* n+N ^ k M n+N j{8M nk 6M* nl ) + M nk M* nl {5M* n+Ntk 8M n+N j) 

k,l=l 



+ M n+N ^M nk (5M* nl 5M* n+Nk ) + M* nl M* n+N ^8M n+N ,5Mnk)\ (Dll) 

In absence of time-reversal symmetry, the contribution from the last two terms inside the bracket is zero (see 
eq.(D7)). However the time-reversal symmetry along with current conservation condition in a region implies M mk — 
^■\m+N\\k+N\- A s a result, the last two terms inside the bracket can be rewritten as follows 



M nk M n+N ,{5M* n+Nk 5M* nl ) = M{ n+NUk+N \M n+N , l [6M n \ k+N \8M* t )8 pi 

M* n+Nik M*JJM^5M^J) = M nt]k+N] M* nl {5M*^ +m 5M nk )5 01 (D12) 
The presence of time-reversal symmetry as well its absence therefore results in 

N 

I 2 - mY ]T S kl [(K +JV , fc M„ w + M nk M* nl ){l + 6 10 )] (D13) 
k,i=i 

The above equation can further be simplified by using cq.(D3) which gives Ajj = 2J2 k =i Mj k M* k — g. Using the 
equality in eq.Q, we get 



\SA n , n+N \ 2 = 8(A nn + A n+N . n+N )5Y (D14) 



1G 



3. Case m = n 



For a clear understanding, here we consider the cases with and without time-reversal symmetry separately. As 
shown below, |<5A„„| 2 turn out to be same in both the cases. 

Region without Time- Reversal Symmetry, (/? = 2): 

As A nn = ELl Efc=l M% k ;a = E r =0,JV ELl EfcLl M lk+r^ we g et 

JV 

,fe+n;si ) (D15) 

n,r 2 si ,S2 fc,J=l 

JV 

,fc+ri;«i ,Z+r 2 ;s 2 ^/c+r! ,Z+r 2 <^Si ,s 2 (D16) 

ri,r 2 si,S2 fc,;=l 
JV 

r s fc=l 

= 8A nn( 5y. (D18) 

Here eq.(D16) is obtained from eq.(D5) by using eq.(24) (or alternatively eq.(D7)). 
Region with Time-Reversal Symmetry 

As M n k-s = (— l) s_1 M| n+ 7v|,|fc+jv|;s under time-reversal symmetry, A nn can be rewritten as A nn = 

Er=0,JV ELl EfcLl M l+rM;s- This S iveS 

JV 

SA nn 6A* nn = 4 E E Mn + r * < k ^ M n+r 2 •*■<» ( SM n+n ,k;«i SM n+r2 ,l;s 2 ) (D19) 

n,r 2 si,s 2 fc,i=l 

JV 

+^ 2 ,i;s 2 

+ (-1)" 2 <5 riir2+ jvc>fc,;+jv<5si,s2] (D20) 

i"i,r 2 si,s 2 Jc,i = l 
JV 

= W EEEt M ^M + (-l) , " 1 ^+ r+ «,Hi;»^+r,fc;»] (D21) 

r s fe=l 
JV 

= 8<5r EEE M «+^ (D22) 

r s fe=l 

= 8A nn 6Y (D23) 

here eq.(D22) is obtained from eq.(D21) by using the relation M n+r+ N.k+N-.s = (— l) s 1 M n+r ^ s - 
The eqs.(D9),*(D14) and (D23) can be combined together to give eq.(25). 



APPENDIX E: PROOF OF EQUATION (10) 

The matrix A is Hcrmitian in nature. Let the inverse pair eigenvalues of A, at value Y of the complexity parameter, 
be given by x n and x" 1 , n = 1 — > N at parametric values Y". A small change SY in parameter Y changes A and 
its eigenvalues. By considering matrix A + 5 A in the diagonal representation of matrix A, a small change 5x n in the 
eigenvalues can be given as 

6x n = 6 A m + E + °((^-) 3 ) ( E1 ) 

where A mn = x n (5 mn at value y of complexity parameter. This further gives, 
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2JV 



Sx n — SA nn + y y 



2N 



-2 1 A nn + 2 ]T (/5 + (4 - /3)<5 

m—l,m^n 



A a- A 



SY 



2N 



-2 1Xn + 2 Y, G 9 + ( 4 - P) 5 ra,n +N ) 

= F n (x)5Y 



x n -\~ x m 



SY 



Here eq.(E3) in the above has been obtained from eq.(E2) by using eq.(25). Similarly, upto first order of SY, 



SXnSXva — 8A nn 8A mm — &A nn Snrn,8Y — QXn&nvn 

SY 

The eq.(E5) can now be used to obtain first moment of eigenvalues A„ = ^ Xn 4x 1 " 1 of the matrix B. 



5\ n 



4r 2 = ; K(A)5F 



(E2) 

(E3) 

(E4) 
(E5) 

(E6) 
(E7) 



where 



F n = 



N 



-2 1Xn + 2/3 

m—l : m^n 



(x n + X m ) (x n X rn ) + {x n Xin ) + X m ) X n + X n 



= ^ [2 7 ,„VA„(1 + A„)- 4(2A„ + 1) - 4/3 ]T Mi+M] 
.x fJ 



-En 



where i?„ is given by eq.(28). The average of the second moment (SX n ) 2 can similarly be calculated, 



SX m S\ 

n 



(xl - 1) (xl 1) 

4x™ 4x„ 



fenfem = 8 — — ^ S nm SY 



(E8) 



(E9) 
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